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We report on the first experimental realization of optimal linear-optical controlled phase gates for arbitrary 
phases. The realized scheme is entirely flexible in that the phase shift can be tuned to any given value. All such 
controlled phase gates are optimal in the sense that they operate at the maximum possible success probabilities 
that are achievable within the framework of any postselected linear-optical implementation. The quantum gate 
is implemented using bulk optical elements and polarization encoding of qubit states. We have experimentally 
explored the remarkable observation that the optimum success probability is not monotone in the phase. 
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Linear-optical architectures belong to the most prominent 
platforms for realizing protocols of quantum information pro- 
cessing lllllj]- In small-scale applications of quantum infor- 
mation, such as in quantum repeaters, they will quite cer- 
tainly play a key role. Unsurprisingly, a significant research 
effort has been dedicated in recent years to experimental re- 
alization of universal linear-optical quantum gates. Linear- 
optical quantum gates are probabilistic by their very nature 
111]. Therefore, the exact trade-offs between properties of a 
gate, such as entangling power, and its probability of success 
are in the focus of attention. 

In this work, we explore this trade-off for the first time 
experimentally. We present data from an experimental real- 
ization of a linear-optical, post-selected controlled phase gate 
implementing the following operation on two qubits: 



|0,0) ^ 2io,o|0,0) = |0,0), 

|0,1) ^ ^io,i|0,l) = |0,l), 

|1,0) ^ ^ii,o|l,0) = |l,0), 

|1,1) ^ ^ii,i|l,l) = e^^|l,l), 



(1) 



for an arbitrary given phase (^ G [0, tt]. It is key to this exper- 
iment that this angle can be chosen in a fully tunable fashion, 
hence adding a flexible scheme to the linear optical quantum- 
information-processing toolbox. 

Controlled phase gates are important members of this tool- 
box. For example, they play a key role in the circuit for 
quantum Fourier transform [3]. They are entangling quan- 
tum gates in general and, together with single-qubit opera- 
tions, they form a universal set for quantum computing. No- 
tice that the controlled-NOT gate can be obtained by applying 
a Hadamard transform to the target qubit before and after a 
controlled phase gate with phase shift tt. 

Previous experimental work was devoted to the linear- 
optical realization of a special case of the controlled phase 
gate with the fixed phase (/:? = tt [4]. Ref. [5] presents an ex- 
periment with phases different from tt, but with a non-optimal 
probability of success. The optimal success probability has 
recently been identified theoretically in Ref. 16]. This opti- 



mum probability we have indeed reached in the experiment 
described in this Letter. We observe the quite remarkable 
trade-off between the entangling power of the gate and its suc- 
cess probability, which is — surprisingly — not monotonous in 
the phase on [0, tt]. The success probability decreases rapidly 
for small phases, but remains almost constant for phases be- 
tween 7r/4 and tt. This experiment is hence expected to be 
both interesting conceptually as well as technologically, since 
a fully tunable bulk linear-optical architecture is presented, 
uplifting tunable schemes for quantum state preparation \m to 
the level of quantum information processing. 

Theoretical framework: We start by introducing the the- 
oretical framework underlying the experiment. For post- 
selected linear optical gates, the beam splitter matrix A de- 
scribing a general linear optics network is constrained by the 
action of the gate ([B as 



perA[ci,Cj|c/e,Q] = UijSi^k^j^i, 

where i^j^k^l = 0,1. The left hand side is the permanent of 
a matrix filled by matrix elements of A [8]. cq = (0, 1), ci = 
(1,0) are vectors describing the usual dual-rail encoding into 
Fock states on two modes [9]. Due to post-selection, only out- 
comes in the computational subspace span{|c/e, ci)} are con- 
sidered, giving rise to 16 quadratic equations in the matrix el- 
ements of A. The solutions are (up to mode-permutations) all 
of the form A = IL2 5 with B being a 2 x 2-matrix describ- 
ing the interaction of the two logical- 1 modes while the oth- 
ers are just by -passed. Using the singular value decomposi- 
tion, this matrix can be expressed sls B = ]/diag(cr+, a-)W, 
where the unitaries V and W have immediate physical mean- 
ing: they describe beam splitters between the two modes. Af- 
ter rescaling the diagonal matrix by the largest singular value 
a-^, an interpretation of diag(l, a-/a^) can be given as well: 
One mode is left unaltered while the second one is damped 
by 6 = g-Jgj^. The success probability is then given by 
Ps = cr^^. It turns out \^ that for Lp G [0,7r] this optimal 
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FIG. 1 : Conceptual scheme of the gate. Vertically (F ) and horizon- 
tally {H) polarized components of the same beam are drawn sep- 
arately for clarity. In polarization beam splitters PBSl and PBS2 
the vertical components are reflected. Half-wave plates HWP6 and 
HWPc act as "beam splitters" for V and H polarization modes. 
Fl and F2 are filters (attenuators), Fl acts on the both polarization 
modes, F2 on the H component only. Phase shifts 0+ and 0- are in- 
troduced by proper path differences in the respective modes. HWPa 
and HWPd just swap vertical and horizontal polarizations. In the 
final setup they are omitted for simplicity and the second qubit is 
encoded inversely with respect to the first qubit. 



success probability takes the form 



Ps(v^) 
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and W 



/V2, V = iy-Miag(e^^+,e*^-), 



-1 -1 
1 -1 

where the phases in the lower Mach-Zehnder interferometer 
(between HWP6 and HWPc, see Fig.[T]) are defined by 



arccot 



4 V 4 



Phase shifts are applied to each arm and one arm is damped 
by an attenuator (neutral-density filter) which is characterized 
by an amplitude transmissivity of 



1 



2sinf 
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1 + 2 sin f + 2(2 - 2 cos (/:)) V4 cos 
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The remaining two attenuators in the upper beams are used 
to damp the amplitude of the by-passed modes to compensate 
for the overall losses in the lower beams. Their amplitude 
transmissivity reads 7 = Ps • 

Details of the experiment: As the starting point of this ex- 
periment we generate a pair of photons in the process of type-I 
spontaneous parametric down-conversion. The laser beam of 
250 mW of cw optical power emitted by Krypton-ion laser at 
413 nm impinges on the LilOs crystal. Pairs of photons at 
826 nm are collected using single mode fibres serving also as 
spatial filters. Subsequently, polarization controllers are em- 
ployed to adjust the horizontal polarization of the photons. 

The half-wave plates (HWP) and quarter-wave plates 
(QWP) in the input arms (see Fig. ^ are used to set the in- 
put states. Subsequently, the photons are superposed on the 
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FIG. 2: Scheme of the actual experimental setup (see text for details). 



first polarizing beam splitter PBSl which transmits horizon- 
tal and reflects vertical polarization. Due to imperfections the 
transmissivity for horizontal polarization is only 95 % (the re- 
maining 5 % are reflected). Polarization beam splitters also 
introduce parasitic phase shifts between vertical and horizon- 
tal polarization components. After leaving the PBSl the pho- 
tons in the upper arm are subjected to the action of half- wave 
plate HWP21. When set to 22.5 deg it performs the transfor- 
mation \H) ^ (\H)^ \V))/V2, \V) ^ (\H)- \V))/V2, 
where \H) and \V) denote horizontal and vertical polariza- 
tion states, respectively. The lower arm is also equipped with 
a half wave-plate (HWPll) but it is set to zero (its presence 
just guarantees the same optical paths, dispersion effects, etc. 
in the both arms). Behind the wave plates there are the beam- 
divider assemblies BDAl and BDA2. They consist of two 
beam dividers (BD) splitting and subsequently rejoining hori- 
zontal and vertical polarizations. BDA2 is equipped with gra- 
dient neutral-density filters F21 and F22 (see Fig. [21). This 
way one can perform arbitrary polarization sensitive losses. 
BDAl is used just to equilibrate the beam position and the op- 
tical length of the both arms. After leaving the beam-divider 
assemblies the photons propagate through half-wave plates 
HWP12 and HWP22. HWP22 is set to 22.5 deg reversing 
thus the transformation imposed by HWP21. HWP 12 is set 
to 45 deg to compensate for the polarization flip between the 
H and V polarizations performed by BDAl. The lower arm 
is equipped with a gradient neutral density filter Fl to apply 
polarization independent losses. The gate operation itself is 
completed by overlapping the photons on the second polariz- 
ing beam splitter PBS2. To be able to perform complete state 
and ihtroby process tomography we employ polarization anal- 
ysis in the both output arms. The analysis consists of QWPs 
and HWPs followed by polarizing beam splitters, cut-off fil- 
ters and single mode fibres leading to single photon detectors. 
Gate operation: The setup is then adjusted to perform the 
gate operation. First we set filters F21 and Fl to introduce 
the required losses. After that the wave plates HWP21 and 
HWP22 are set to 22.5 deg. The phase in the beam divider 
assembly BDA2 is set to maximize the visibility of the inter- 
ferometer formed by PBSl and PBS2. The precise tuning of 




O.OStt 
0.1257r 
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0.94 
0.95 
0.91 
0.84 
0.86 
0.84 
0.84 



0.96 
0.96 
0.90 
0.88 
0.89 
0.87 
0.86 



0.84 
0.91 
0.77 
0.73 
0.81 
0.63 
0.71 



0.96 
0.96 
0.95 
0.90 
0.90 
0.90 
0.92 



0.87 
0.87 
0.87 
0.67 
0.76 
0.71 
0.83 



Ps,obs 

0.87 ± 0.08 
0.37 ± 0.05 
0.19 ±0.03 
0.11 ±0.02 
0.09 ± 0.01 
0.08 ± 0.01 
0.12 ±0.01 



Ps,th 

1.00 
0.36 
0.21 
0.13 
0.09 
0.09 
0.11 



TABLE I: Process fidelities {F^), average (Fav) and minimal (Fmin) 
output-state fidelities, average (T^av) and minimal (T^min) output- 
state purities and actually observed (ps,obs) and theoretically pre- 
dicted (ps,th) success probabilities for different phases {Lp). 





FIG. 3: Matrix representation of the CP map characterizing the oper- 
ation of the controlled-phase gate with Lp — n. The left panel shows 
the real part of the reconstructed process matrix, the right one dis- 
plays the real part of the ideal theoretical CP map. Imaginary parts 
are negligible (zero in ideal case). The process fidelity F^ = 84%. 



the gate is then performed by switching between the inputs 
\Hi^R2) and | Fi , i?2), where indices 1 and 2 denote the input 
modes and R stands for the right circular polarization. Using 
the circular detection basis in the second output arm we can 
observe the phase applied by the gate when the polarization 
of the first input photon flips from \H) to \V). In this config- 
uration we also tune the phase shift inside the beam divider 
assembly BDA2 and the phase shift between the two arms of 
the Mach-Zehnder interferometer formed by PBS 1 and PBS2. 

Results: Gradually we have adjusted the gate to apply 7 
phases in the range between and tt. Each time we have per- 
formed complete process tomography and estimated the pro- 
cess matrix using the maximum likelihood method. Eidelities 
of the process lie in the range from 84% to 95% (see Tab.!]). 
Figs. [3] and [4] show examples of experimentally obtained pro- 
cess matrices and their theoretical counterparts for (/:) = tt and 
7r/2. 

Eor each selected phase we simultaneously measured two- 
photon coincidence counts between detectors D1H&D2H, 
D1V&D2V, D1H&D2V, and D1V&D2H, each for 3 x 3 com- 
binations of polarization measurement bases in the output 
arms. This amounts to measuring projections onto horizon- 
tal/vertical, diagonal/anti-diagonal and right/left circular po- 
larizations. The diagonal and anti-diagonal linear polarization 
states are defined as I D) = {\H)^V))/^/2^n&\A) = (\H)- 
\V))/y/2 and the right- and left-handed circular polarization 



FIG. 4: Choi matrices for the gate with Lp — 7r/2. The left top 
panel shows the real part of the reconstructed process matrix while 
the left bottom one displays its imaginary part. The process fidelity 
F^ — 86%. The two right panels show the ideal matrix. 



states read I i?) = {\H) ^ i\V)) / ^/2, \L) = {\H) -i\V))/V2. 
The unequal detector efficiencies were compensated by proper 
re-scaling of the measured coincidence counts |K)I]. Each 
measurement was done for 36 different input product states. 
Namely, for 6 x 6 combinations of polarization state vec- 
tors \H), \V), \D), \A), \R), and \L) of each input photon. 
This complex measurement provided us with tomographic ally 
complete data enabling us to fully characterize the imple- 
mented operation by quantum process tomography ||jjj, \l^ 
as well as to reconstruct density matrices of output states for 
each used input state. 

Active stabilization: Each setting of an input and output 
polarization basis was preceded by an active stabilization. For 
the purpose of the stabilization the fixed input state and output 
detection basis were always used. In this setting the visibility 
in the interferometer formed by PBSl and PBS 2 was mea- 
sured. If this visibility was lower than a selected threshold 
then the positions of MTl (interferometer lengths) and MT2 
(dip position) were optimized and the phase drift was com- 
pensated. Finally the required polarizations were set and data 
were accumulated within 5 s. 

Process tomography: Any quantum operation can be fully 
described by a completely positive map and — according to the 
Jamiolkowski-Choi isomorphism — represented by a positive- 
semidefinite operator y on the tensor product of input and out- 
put Hilbert spaces U3 . In our case xisal6xl6 square 
matrix. From the measured data we can reconstruct v for any 
setting of cp using maximum likelihood estimation |[l3,[lJ] . To 
quantify the quality of the operation we calculate the process 
fidelity, if Xid is a one-dimensional projector its common def- 
inition is F^ = Tr[xXid]/(Tr[x]Tr[xid]). Here Xid represents 
the ideal transformation corresponding to the controlled-phase 
gate. Specifically, 

Xid= J2 \i,J){k,l\^U\i,j){k,l\U^, 
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FIG. 5: Success probability of the gate. 



where U stands for the unitary operator on two qubits defined 
by Eq. ([T]). We have also reconstructed the density matrices 
of output two-photon states corresponding to all product in- 
puts |j, /c), j, /c G {H^V^D^A^R^L}. This was done for 
all values oi (p. An important parameter characterizing the 
gate performance is the fidelity of output states pout defined 
as F = (V^outlPoutlV^out), where |V^out) = ^IV^in) and |V^in) 
is the input state vector. Table U contains the average and min- 
imal values of state fidelities for different phases. Fidelities 
Fav are averaged over all output states corresponding to our 
36 input states; Fmin denote minimal values. Another im- 
portant characteristics is the purity of the output state pout, 
defined as 7^ = Tr[pQ^^]. If the input state is pure the output 
state is expected to be pure as well. The average and minimal 
purities of output states are also given in Table H 

Trade off in success probabilities: The most impor- 
tant result of this paper — aside from the technological 
implications — is the experimental verification of the trade-off 
between the phase shift applied by the gate and the corre- 
sponding success probability of the gate. We have estimated 
the success probability for each value of the selected phase 
shifts. It was calculated as a ratio of the number of successful 
gate operations per time interval and the number of reference 
counts during the same interval (measured with no filters and 
with the wave plates set to 0). We have determined the suc- 
cess probability for the all selected input states. These proba- 
bilities were averaged and the standard deviations were calcu- 
lated. Notice that the calibration measurements collect coin- 
cidence counts behind the setup (using the same detectors as 
in the subsequent measurements), thus all the "technological" 
losses in the setup and low detector efficiencies are included in 
the calibration. Therefore the estimated success probabilities 
are not burdened by these "technological" losses. They can 
be compared with the theoretical predictions in Tab. U and in 
Fig.O One can see a very good agreement with the theoretical 
prediction. 

Conclusions: We have built the first implementation of the 
tunable linear-optical controlled phase gate which is optimal 
for any value of the phase shift. Changing the parameters of 



the setup the gate can apply any phase shift from the interval 
[0, tt] on the controlled qubit. We have thoroughly tested the 
performance of the gate using full quantum process tomog- 
raphy. Obtained process fidelities range from 84% to 95%. 
We have also experimentally verified that all our controlled 
phase gates are optimal in the sense that they operate at the 
maximum possible success probabilities that are achievable 
by linear-optical setups. The experimental verification of this 
trade-off between the phase shift applied by the gate and the 
corresponding success probability of the gate is the most no- 
table result of our work. It demonstrates the contra-intuitive 
fact that the optimal success probability is not monotonous 
with the phase shift increasing from to tt. It is the hope 
that the flexible tool established here proves useful in devis- 
ing further optical linear optical circuits for optical quantum 
information processing and that ideas developed in this work 
find their way to realization in fully integrated optical archi- 
tectures. 
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